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ABSTRACT 

We present a multi-wavelength photometric and spectroscopic study of a 
newly discovered candidate cluster [DBS2003] 45. Our H, Ks photometry con- 
firms that [DBS2003] 45 is a cluster. An average visual extinction Ay~7.1±0.5 is 
needed to fit the cluster sequence with a model isochrone. Low resolution spec- 
troscopy indicates that half a dozen early B and at least one late O type giant 
stars are present in the cluster. We estimate the age of the cluster to be between 
5 and 8 Myr based on spectroscopic analysis. Assuming an age of 6 Myr, we fit 
the observed mass function with a power law, N(M) oc M _r , and find an index 
T ~ 1.27±0.15, which is consistent with the Salpeter value. We estimate the total 
cluster mass is around 10 3 M Q by integrating the derived mass function between 
0.5 and 45 M . Both mid-infrared and radio wavelength observations show that 
a bubble filled with ionized gas is associated with the cluster. The total ionizing 
photon flux estimated from radio continuum measurements is consistent with the 
number of hot stars we detected. Infrared bright point sources along the rim of 
the bubble suggest that there is triggered star formation at the periphery of the 
HII region. 

Subject headings: Stars: evolution, Stars: luminosity function, Galaxy: open 
clusters and associations 
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Introduction 



Stellar clusters serve as ideal laboratories to study the evolution of stars of different 
masses, as the stars in a cluster are formed at the same time from the same molecular cloud. 
This assumption predicts that the stars are at the same distance, have practically the same 
initial chemical abundances, and evolve in a similar physical environment. Such a coeval 
population can be used to make critical tests of theoretical stellar evolution models. Young 
clusters are of particular interest for this purpose because most of their massive stellar 
members are still in main sequence (MS) or post-MS stages before supernova explosions. 
Studies of such stars can improve our understanding about the formation and evolution of 
massive stars, as well as the evolution of whole clusters in general. 

Stars generally form with a frequ ency that decreases with increasing mass for M>1 
M Q , i.e. r=d(logN)/d(logm) ~ -1.35 JSalpeterlll955l : IScaldfl998l : iKroupafeoOll l2002h . This 
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(2002 


) and 


Kim et al. 



( 120061 ) claim a significantly shallower slope for the Arches cluster. This result, for one of the 
most massive and the densest cluster in the Galaxy, suggests that more massive and dense 
clusters have a shallow IMF. Alternatively, perhaps the Galactic center f avors sha l low IM Fs 
( lMorriall993l ) . Since the existing IMF slope of -1.35 was first measured by lSalpeterl (119551 ) for 
field stars with masses lower than 10 M and field stars have a MF which is different from the 



IMF of clusters, the claim of a universal IMF seems questionable. IZinnecker fc Yorkd (120071 ) 
commented that, although massive stars can form locally in the field, evidence suggested 
that massive stars form preferentially in dense clusters or OB associations. Additionally, 
mass segregation can cause massive stars to sink toward the cluster center and the cluster 
continuously loses low mass stars to the field. These effects would modify the cluster MF so 
that it becomes shallower as the age of the cluster increases. The ultimate resolution to this 
question will come from measurements of the IMF in stellar clusters throughout the Galaxy. 
The sample of newly- identified clusters provides a fresh opportunity for this purpose. 

There are le ss than 1600 open clusters discovered at optical wavelengths in our galaxy 
( IDias et al.ll2002l ). It is reasonable to believe that a large number of such clusters are missed 
by optical observations due to substantial line-of-sight interstellar extinction in the Galactic 
plane. The advance of infrared technology allows surveys to be carried out in near infrared 
(NIR) and mid infrared (mid-IR) wavelengths to circumvent the problem of high extinction. 
Systematic searches for Galactic open cluster candidates hav e been carried out in infrared 



bands and more than 1000 new candidates were discovered ( Bica et al.ll2003l ; iDutra et al. 



2003al lbl: 1 Froebrich et al.l 120071 ). although many of them remain unconfirmed and extensive 
efforts in follow-up observations are needed to determine their true identities.. 
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Here we present our photometric and spectroscopic analysis of [DBS200 3] 45, a candidate 
cluster chosen from such an infrared cluster catalog (IDutra et al.l l2003al ). This target is 
located at (1, b) = (283. °88, -0.°91) and is about 3.5 degrees on the sky to the west of 
the Carina nebula (NGC3372), suggesting that the candidate cluster is likely located in the 
Sagittarius-Carina spiral arm. In § 2, we describe the origin of data that are used in this 
work. Our data reduction procedure is described in § 3. In § 4, we present the results of our 
analysis. Some issues related to our observation are discussed in § 5. Finally, we summarize 
our findings in § 6. 



2. Observations 

Photometric observations of [DBS2003] 45 were carried out on the ESO 3.8m New 
Technology Telescope (NTT) with SofI in 2008 February (ESO proposal ID 080.D-0470 ). 
SofI is an infrared spectrograph and imaging came ra equipped with a Ro ckwell Hg:Cd:Te 



1024 x 1024 HAWAII array of 18.5 micron pixel-size (IMoorwood et al.lll998l ). It takes images 
with a plate scale of 0.288"/pixel. H and Ks band images of the target field were obtained 
in a dither pattern in which candidate clusters were shifted to one of four corners of the 
detector each time. We used discrete integration times (DITs) of 2s, coadded 10 times, in 
4 dithers. Multiple images were taken at each dither location, resulting in a total 160s of 
integration time for the location of the cluster at each waveband. 

We also obtained low-resolution spectroscopy of several stars in the field of the cluster. 
We used the 'Red' grism in conjunction with the 0.6" slit, giving a spectral resolution of 
R ~ 900 in the wavelength range 1.9 - 2.4//m. We employed DITs of 40-60s, repeated 
twice at each position. After each integration the target was nodded 20" along the slit. 
Each object was observed in 8 nod positions, hence total integration times were 16mins per 
object. For fainter objects (Ks > 10) we increased the number of nod positions to 16 in 
order to improve signal-to-noise. To correct for atmospheric absorption, the solar-type star 
Hip059532 was observed immediately after observing the cluster stars. Dark-frames and flat- 
field observations were made during twilight. For wavelength calibration purposes a neon 
arc lamp was observed. 

Additionally, images of the field were taken in four bands at 3.6, 4.5, 5.8, and 8 //m on 
February 2007 using the infrared camera - IRAC on the Spitzer space telescope (Program 
ID: 30734). During the observations, we used the 5 position Gaussian dithering pattern. 
The resulting images of each band cover an area of 5' by 5' with a pixel size of 1.2" in both 
directions of RA and DEC. The total integration time on target is 471s. 
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In order to study the nebulosity surrounding the cluster and its association with the clus- 
ter, we use H Q , mid- infrared maps, and radio survey d ata. H a line observations were taken 
from superCOSMOS survey archive (IParker et al.ll2005f). Radio flux is ob tained from Parkes- 
MIT-NRAO 4850 MHz southern sky survey catalog (I Wright et al.lll994f ). We also obtained 
from the archive an 8 /im image of a large area including our field taken by Spitzer/IRAC 
(Program ID: 40791). The mid-IR image and the H a emission map will be compared in 
order to address the association between the cluster and the bubble. 



3. Data Reduction and Analysis 

Our photometric and spectroscopic data sets come from different sources. Therefore, 
the reduction procedure is slightly different for individual data set. The SofI data set in- 
cludes both photometric and spectroscopic observations. Here we describe their reduction 
procedures separately. 



3.1. SofI Data 

3.1.1. Photometry 

Our Soft photometric data consists of two bands, H and Ks, with eight partially over- 
lapping images in each band. The flat-fielding is done using facility provided dark and flat 
frames. A group of 8 images that were obtained closest in time are scaled according to the 
integration time and median-combined to create a sky background image. This background 
image is scaled and subtracted from each flat-fielded image of the target. The resulting im- 
ages are re-sampled, cross-correlated, shifted, and coadded to create a mosaic. We register 
resulting mosaics with downloaded 2MASS mosaics of the same field in the same photomet- 
rical bands to correct for astrometric errors in our observations. We manually pick a group of 
relatively bright and isolated stars that can be seen in both SofI and 2MASS mosaics. Using 
the 2MASS coordinates of these bright stars, we compute the geometric transformation co- 
efficients between the coordinates of observed images and the 2MASS images. The observed 
images are re-mapped onto the 2MASS grid to produce the final images for photometric 
analysis. The final Ks band image is shown in Figure [TJ 
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3.1.2. Calibration 

The aperture photometry of point sources are obtained for the SofI data using an IDL 
adapted version of the DAOPHOT. Since all images are re-mapped onto the 2MASS grid, 
the magnitudes of different bands can be associated using stellar locations. The same group 
of stars that were used to register the SofI mosaics are used to calibrate detected H and 
Ks band magnitudes. Figure EJ^left) shows corresponding instrumental and 2MASS magni- 
tudes of these stars with plus signs and diamonds representing H and Ks band magnitudes, 
respectively. It clearly shows that linear relations (solid lines) exist between two sets of 
measurements. 

To account for any color term effect present in the data, the following linear bivariate 
model between the instrumental magnitudes and the 2MASS magnitudes is assumed: 

h = a h + b h H + c h Ks (1) 
ks = a k + b k H + c k Ks (2) 

here, h and ks are instrumental magnitudes of the picked stars. H and Ks are the cor- 
responding 2MASS magnitudes. The model coefficients are determined using a regression 
method. Once all coefficients have been determined, the following equations are used to 
compute calibrated magnitudes of all stars in the field 

c k h - c h ks - a h c k + a k c h 

H = u 1 \y) 

T s b k h - b h ks - a h b k + a k b h 

Ks = (4) 

c/A - b h c k 

or 

H = A H + B H h + C H ks (5) 
Ks = A K + B K h + C K ks (6) 

here, coefficients in the capital form (A, B, C) are the results of combining the corresponding 
terms in the equation [3j Table [1] shows the coefficients found by our regression procedure 
and indicates small amount of color effect. Figure [2] (right) shows the corrections between in- 
strumental and calibrated magnitudes. Again, diamonds and plus signs in the figure indicate 
Ks and H band magnitudes, respectively. 



3.1.3. Completeness Test 



The ability to detect stars is impaired when dealing with crowded fields. Faint stars 
are likely missed by a star finding procedure due to the overlapping of stellar point spread 
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functions (PSFs). A completeness test is carried out to correct for such an error and to 
establish the actual magnitude distribution. We add artificial stars to the mosaics and 
apply the same star extracting procedure to see whether or not they can be recovered. Stars 
detected in each mosaic are ranked based on their distances from neighboring stars. The first 
20 unsaturated stars with the largest separations are picked. Each star's PSF is extracted 
from the mosaic and the local sky background is subtracted. Finally, the individual PSFs 
are re-sampled, shifted, and added to create a template PSF. 

We divide the magnitude range from 9 to 19 magnitudes into bins with a binsize of 
1 magnitude. Within each bin, we create 50 fake stars with their magnitudes uniformly 
distributed across the bin. The images of the stars are created by scaling the template PSF 
and added to the input mosaic at random locations within 43" from the center of the cluster 
field. The resulting mosaic is then treated with DAOPHOT in the exactly same manner 
that is used to find stars in the original extraction. If a star is found at the location where 
a fake star is added and the difference between the input and output magnitudes is less 
than 0.5 mag, we define that the fake star is recovered. Otherwise, we define that the added 
artificial star is missing. We repeat this experiment 10 times and generate 500 artificial stars 
total which are evenly distributed over the magnitude bins. Figure [3] shows the curve of 
completeness from our simulation, which suggests that we achieve the 50% completeness at 
Mk s ~17 (mag) at the center of the field. 



3.1.4- Spectroscopy 

We first subtracted dark frame from each science frame to minimize fixed-pattern noise 
of the detector. Observations at two complimentary nod-positions were then subtracted from 
each other to remove atmospheric OH emission lines. Repeating observations were averaged 
to reduce noise. Afterwards, spectra were extracted from the average frame by summing 
the rows across the width of the spectral traces. A spectrum of the arc lamp at the same 
position on the chip was also extracted using the same extraction procedure. Wavelength 
calibration was achieved by fitting the locations of the arc lines with a 4th-degree polynomial. 
Atmospheric absorption in the stellar spectra was removed using the telluric standard star, 
which had first been divided through by a synthetic spectrum of a GV star. Finally, stellar 
spectra were normalized by dividing through by their median continuum values. The reduced 
spectra are shown in Figure HI The coordinates of the stars are listed in Table HJ From flat 
regions of continuum, we estimate the signal-to-noise of the spectra are 80 or better. 
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3.2. Spitzer IRAC Photometry 

The basic calibrated data (bed) was downloaded and used to create new mosaics using 
Spitzer's Mopex tool before the mosaics were treated with DAOPHOT. Point spread function 
(psf) photometry was performed in order to achieve deeper photometric detection. After 
every round of star finding process, the detected stars were subtracted from the mosaic 
to create a residual mosaic. Additional rounds of star finding were then performed on the 
resulting mosaics until no additional stars were found in the final image. Finally, all detected 
point sources were collected and cross-correlated with 2MASS point sources for completeness. 

4. Results 

4.1. Color Magnitude Diagram 

Figure EKright) shows color magnitude diagrams (CMDs) of the candidate cluster before 
(left panel) and after (middle panel) field contamination has been removed and the CMD of 
a chosen control field (right panel) for our SofI data. The image on the left demonstrates 
how the regions for the cluster and the control field are defined. A circular region with a 
radius ~0.72 arcmin is chosen for the cluster and it includes the densest part of the observed 
field. An annulus region outside the circular region is chosen as the control field, which has 
an outer radius about 2.4 arcmin. Both regions are centered at RA(2000)=10 /l 19 m 10.5 s , 
Dec(2000)=-58°02'22.6". The CMD of the control field is scaled based on the area ratio of 
regions for the cluster and the field and is subtracted from the CMD of the cluster+field. 
The resulting decontaminated cluster CMD is shown in the middle panel of Figure fright), 
which demonstrates clearly a cluster sequence with majority of stars along a vertical line with 
H-Ks~0.45. Among 175 remaining stars, the brightest one has Ks~10. Since a few stars are 
located above our saturation limit (dashed lines in the CMDs) , we have replaced presumably 
saturated magnitudes with the corresponding values in the 2MASS catalog, which does not 
significantly change the locations of those stars on our CMDs. 

4.2. Spectroscopy 

Low-resolution spectra of a few relatively bright stars present in the field of the candidate 
cluster are shown in the right panel of Figure H] with their locations indicated in a Ks band 
image of the field on the left together with other stars remaining in the cluster CMD. The 
spectrum of the brightest star in the field #8 shows CO bandhead absorption. The CO 
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equivalent width is consistent with the star being either a late-type giant (M5-7III) or a 
supergiant of an earlier type (K2-4 I) . The spectrum of the star #6 shows faint emission at 
the wavelengths of the Br7 line and the Hel 2.06 /mi line and absorption in the Hel 2.11 /mi 
line. All other spectra show weak absorption at the wavelengths of the Br7 line, the 2.06 /mi 
and 2.11 /mi Hel lines. None of these spectra shows significant contribution from the 2.19 /mi 
Hell line. 



According to the atlas provided by lHanson et al.l (120051 ). stars earlier than 09.5 usually 



show Hell 2.19 /mi line in absorption in their spectra. This feature should be as strong 
as the Br7 line absorption in the early O type stars, fades away after 07 while Br7 gets 
stronger. The spectra of late O and early B type stars should also possess the Hel 2.06 /tm 
or 2.11 /mi line in absorption. These qualitative comparisons suggest the hot stars shown in 
Figure H] have spectral types later than 08, probably early B types. Particularly, star ^6 has 
a spectral type earlier than BO. The presence of the Br7 line emission indicates a spectral 
type around 09. Accurate spectral types need to be determined by future high resolution 
spectroscopic observations for these stars. 



4.3. Extinction 



The intrinsic H-Ks colors of most main sequence and evolved stars are approximately 
zero. This allows us to estimate the average foreground extinction from the observed H-Ks 
colors by assuming that the observed H-Ks values of stars in a cluster sequence are entirely 
due to the foreground extinction. In this way, we derive a foreground visual extinction 



Ay~7.1±0.5 assuming a standard extinction law (IRieke fc Lebofskyl 119851 ). 



Additionally, patchy extinction may be present at the observed field because the distri- 
bution of stars is not symmetric. A hole of stars can be seen at the west side of the identified 
cluster and more stars seem to concentrate toward the southeast side of the cluster. 



4.4. Distance 



Without high resolution spectra, we have to rel y on measuremen t s in t he literature 
for the kinematic distance to the candidate cluster. iDeharveng et al.l (120051 ) determined 
the kinematic distance of [DBS2003] 45 to be ~4.5 kpc from the Sun based on a radial 
velocity measurement of a molecular cloud-HII region com plex G283.9 -0.9, which seems 
associated with the target based on their loc a tions on sky (lRusseilll2003l ). This is the only 
distance measurement available. iBlum et al.l (120001 ) gave the absolute K band magnitudes 
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for main sequence stars of different spectral types. At a distance of 4.5 kpc, late and 
early B dwarfs will have a in the range of 11 and 14 mags with an assumed foreground 
extinction Ay~ 7.1. The median of these values is at least one magnitude fainter than our 
measurement of ~10.5 mags for those hot stars. This fact suggests that the hot stars in our 
observations are giants or supergiants, which is consistent with the presence of Br7 emission. 

The adopted distance puts the cluster candidate and the Sun at the same radial distance 
from the Galactic center. Therefore, it is reasonable to assume that the members of the 
cluster have similar initial chemical abundances to solar values. For this reason we use solar 
metallicity isochrones in the following analysis. 



4.5. Cluster Age 

In general, it is hard to determine the age of the cluster with only NIR photometric 
data because of the degeneracy of isochrones at the NIR wavelengths. However, the absence 
of early O stars and the presence of several late O or early B type blue giants/supergiants 
in our cluster candidate put s ome constraints on the clu ster age. Based on solar abundance 



Geneva evolutionary models (ILejeune fc Schaererll200ll ). we can determine that the cluster 



must be older than 3 Myr to allow some massive stars to evolve off the main sequence. Our 
isochrone fitting using solar abundance Geneva tracks suggests that for any star earlier than 
07.5 to be in our observed magnitude ranges, the cluster needs to be younger than 6 Myr. 
We applied Ay =7.1 and a distance of 4.5 kpc as suggested in the previous sections during 
the isochrone fitting. Star ^8 in Figure H] also provides some constraint on the cluster age. 
Our spectroscopy suggests that it is either a M5-7 giant or a K2-4 supergiant. The star is 
saturated in our SofI photometry. The 2MASS photometry indicates that the magnitudes 
of the star are H~6.57 Mag and Ks~5.85 Mag. If we take this red giant/supergiant as a 
cluster member, the age of the cluster should be older than 8 Million years. 

To better constrain the age of the cluster, we carry out a Monte-Carlo simulation to 
estimate the number of massive stars that should be observed within our observing range 
(M^<18 and 0.4<H-Ks<0.8) for a cluster. We create a stellar cluster which has an IMF fol- 
lowing the 'Salpeter' law and has a total mass of 10 6 M . In order to improve number statis- 
tics, we use the large cluster mass in the simulation. Afterwards, we scale the resulting num- 
ber down according to an assumed cluster mass. The colors of stars in the cl uster are com 



puted for ages between 1 Myr and 10 Myr according to the Geneva isochrones (jSchaller et al. 



19921 ) . The same distance and average reddening as that derived for [DBS2003] 45 are as- 
sumed for the model cluster. Finally, we count the number o f stars with masse s larger than 



15 M , which is approximately the initial mass of a B0.5 star farum et al.ll200oh . The results 
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are shown in FigureEJ The curves are for a cluster of mass 10 4 M , 2.0 xlO 3 M Q , and 10 3 M . 
The simulation predicts that a cluster with Mciuster ~1O 3 M will have only 2-3 massive stars 
present in our observing range when the cluster is older than 8 Myr. This value is about two 
times smaller than our detection of 7-8 massive stars. This suggests that either the cluster 
mass is over 2000 M Q or the cluster age is below 8 Myr. Using a 9.0 Myr isochrone, we 
estimate that the mass of the cluster is around 1400 M . An alternative 6.0 Myr isochrone 
would suggest a cluster mass of 1170 M Q . Therefore, the choice of the isochrone does not 
cause a significant change in the resulting cluster mass to be above 2000 M . We take this 
as the evidence to support a younger age for the cluster. 

From stellar evolution models, we know that all stars with a mass M<120 M of a 3 Myr 
old cluster should be still on the main sequence. Therefore, there should be a limited number 
of giants or supergiants in the cluster. For stars with an intermediate-mass (M~25 M ) to 
evolve off the main sequence, the cluster should be at least 5 Myr old. Our detection of a 
few evolved stars suggests that [DBS2003] 45 is older than 5 Myr. Therefore stellar models 
support that [DBS2003] 45 is between 5 and 8 Myr old. This analysis also suggests that the 
cool supergiant is not a member of the cluster. 

In Figured we plot a 6.0 Myr isochrone on the top of the CMD of the observed cluster. 
As in Figure [5] and HI the stars with available spectra are indicated with diamonds and 
numbers. Open squares on the top of some stars indicate that IRAC four color photometric 
measurements are also available. In the figure, we also indicate the locations of a few stars 
with different initial masses on the isochrone and the direction of the reddening vector. 

4.6. Mass Function 

The cluster mass function is based on the stellar magnitudes and the chosen isochrones. 
The initial mass of a star can be immediately determined if the star is close enough (within 
5-er uncertainties) to the isochrone in the CMD. It is simply the corresponding initial mass 
of the star with the same magnitude on the isochrone. A few stars are on the left side of the 
isochrone. They could be due to low number statistics. We leave them out in our discussion. 
The initial masses of stars on the right side of the isochrone are found by extrapolating the 
reddening vector back toward the isochrone. Figure [8] (left) shows the mass distribution 
for the final 164 stars with reasonable colors. The Geneva 6 Myr isochrone with solar 
abundances is used to transform colors to masses. The dashed line histogram shows the 
distribution before the incompleteness correction and the solid line shows the result after 
the correction. A set of mass bins equally spaced at the logarithmic scale are chosen to make 
the histogram. A theoretical (Salpeter) mass function of jf- = (jf)~ r with iV(2OM ) = 1 
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and T=1.35 is also indicated with a straight line in the figure. A regression procedure is 
used to fit the power law. The regression only uses the bins with at least 3 counts and with 
a mass bigger than the 50% completeness limit. A slope of r=1.27±0.15, which is close 
to the Salpeter value, is found. We notice that the change of isochrone does not produce 
significant difference in the slope. We estimate the total mass of the cluster by taking the 
summation of the mass function between the mass range from 0.5 M to 45.0 M , which 
results in ~1092 M , ~1170 M , and ~1393 M for the 3 Myr, 6 Myr, and 9 Myr isochrones, 
respectively. 

Our field contamination removed CMD of [DBS2003] 45 has only <200 stars. The small 
size of the sample can result in large statistical uncertainty, especially at the high mass 
end. There are only a limited number of stars in the bins with the largest central masses. 
We investigate the effect of the bin number on the resulting slope. By choose different bin 
numbers, we divide the entire logarithmic mass range into bins of equal sizes. We then 
populate the bins and fit the resulting mass function to a power law. Figure [8] (right) shows 
the result of this analysis. The horizontal axis is the number of mass bins that are used. We 
notice that there is no substantial change in the slope due to the change of the bin size and 
the slope tends to stabilize at ~-1.25 when the bin number approaches 16 from both ends. 
Based on this analysis, we conclude that the slope of the mass function is not very different 
from the Salpeter value. 

Compared to Arches cluster at the Galactic center, [DBS2003] 45 has much lower mass. 
Therefore, our findings can not be used as the evidence to support the statement that the 
mass function of massive clusters is same as less massive clusters or field stars. Nevertheless, 
the current data implies that intermediate or low mass clusters may have a similar mass 
function as field stars. And, the evolutionary effect such as evaporating low mass cluster 
members to the field during the cluster lifetime seems not ch ange th e clust er initial mass 
function significantly. Our result also does not contradict with Morris! ( 19931 ). who suggests 



that clusters near the Galactic center may have shallower IMFs, since [DBS2003] 45 is at a 
similar distance from the Galactic center as the sun. 

However, our finding is based on only one cluster. Low number statistics can be im- 
portant in this case. Observations and analysis of a large sample of clusters are needed to 
answer questions such as whether or not the slope of the initial mass function is a function 
of the cluster mass, the dynamical time and the environment. 
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4.7. Mid IR Excess 

As we mentioned earlier, a certain degree of asymmetry in spatial distribution is observed 
for cluster members. This could suggest patchy extinction toward the field. To confirm this, 
we plot the photometry of stars obtained by Spitzer/IRAC in Figure [9j We can divide 
stars on the IRAC M 3 . 6 _ 4i5 versus M 3 .6-5.8 color-color image into three groups, which are 
indicated by solid-lined boxes. Majority of stars in the lower-left box have M3.6-4.5~O.CI and 
M3.6-5.8~O.Cl (group I). In the second box, stars have almost constant M 3 . 6-4.5 but varying 
M36-58 (group II). In the third box, M36-45 seems to vary linearly with M36-58 and with 
ft^~0.3 (group III). 

Stars in group I include almost all stars on the cluster sequence, particularly those close 
to the assumed isochrone. It is consistent with their small color differences. In the figure, we 
plot the reddening vector with Ax s =6.0 mags assuming the intrinsic M3.6-4.5,o=M3.6-5.8,o=0.0 
and a standard extinction law (A^ocA -1,53 ). With A^ s ~0.112Ay~0.8 mag, the values of 
M3.6-4.5 and M 3 6 _ 5 8 would be small and within the box we plotted. Group II represents 
stars with 5.8 /mi excess. The 5.8 /mi and 8.0 /mi channels of IRAC were designed to 
observe emission from Polycyclic Aromatic Hydrocarbon (PAH) particles in warm interstellar 
medium. The detection of 5.8 /mi excess may suggest a group of stars with strong PAH 
emission from their circumstellar envelopes. One interesting fact about these stars is that 
none of them is present in our cluster CMD. We plot their locations on the top of our Ks 
band image of the field to check the integrity of our results (right panel of Figure [9]), which 
suggests that those 5.8 /mi bright stars (squares) are all at large distances from the cluster 
center, including some areas which are not covered by our NIR observations. Thus, it is 
not a surprise that they are excluded from the cluster population. These stars could be red 
giant stars in the field, which have lower than main sequence stellar effective temperatures. 
Stars in group III distribute approximately along the reddening vector and seems to suggest 
that their colors are strongly affected by varying foreground extinction. Their locations in 
the Ks band map indicate they do distribute unevenly in the cluster, supporting a patchy 
extinction interpretation. 

Stars left unexplained are those faint in the Ks band and heavily reddened in Figure [7J 
None of them have an IRAC correspondent. They are probably too faint for IRAC to 
pick up since the detection limit of IRAC photometry is Ks~15.3 mags based on our cross- 
correlation result between the IRAC catalog and the SofI catalog. Their uneven spatial 
distribution in the cluster demonstrated by Figure H] is similar to that of stars in group III 
in the Figure [91 Therefore, the extra red colors of these faint stars may also be caused by 
foreground extinction, although we could not rule out the possibility that they have intrinsic 
colors similar to those in group II. Deeper photometry is needed to clarify on this. 



13 



Discussion 



5.1. The Cluster's Ionized Bubble 



Ionized gas has been observed around massive clusters. In some cases, well-defined 
bubbles of ionized ma terial and warm dust are seen in radio and infrared wavelengths 
( IDeharveng et al.ll2005l ). It is suggested that dusty material around clusters is heated by 
energetic photons from cluster members to high temperatures and emits at infrared wave- 
lengths. By examining the properties of radio/IR continuum from the heated material, one 
can check the consistency between the total energetics and the assumption about the cluster 
members. 

We over-plot Spitzer IRAC 8.0 fim emission contours on the top of superCOSMOS H a 
emission grey scale image (Figure [TD], left). A well-defined cavity of 8.0 /im emission can 
be seen around relatively concentrated H a emission. The location of the identified cluster 
inside the bubble suggests that the cluster is the source of the ionizing photons. Radio VLA 
4850 MHz continuum ob servations showed t hat the total continuum flux density of the HII 
region is 4519±99 mJy (jWright et al.lll994j ). In a thermal situation, radio continuum flux 
luminosity (Lf/) and flux density (Fj /) are related to the total Lyman continuum photon 
flux, Q Ly ( Rvbicki fc Lightman 19861 ): 



Lf/ 



Qhy 



4nD 2 F/, f 



Q.8xlO~ 38 T°- 5 e~ h ^ kT g ff 



N P K-dV 



(7) 



a B / N e NdV. 



Here D is the heliocentric distance to the HII region, T~10 4 K electron temperature, g//=5.1 
the Gaunt factor, N e and N{ the electron density and the proton density, and as the re- 
combination coefficient. The integral is over the entire volume of the ionized region. We 
can derive the ionizing photon luminosity to be Ql,,~8.2x10 48 sec -1 . This luminosity cor- 
responds to the total ionizing luminosity of one 08.5 giant star or four BO type giants 
(ISchaerer fc de Koterlll997l ; ISternberg et al.l 120031 ) . This is consistent with our detection of 
~6 early B giants plus one 09 giant. The total e nergy budget is also consistent with the 
value from far-IR observations. IDeharveng et al.l (120051 ) reported that the far-IR fluxes of 
IRAS point source associated with the dust ring indicates that the total energy required to 
power the IRAS source is equivalent to the total luminosity of 10 B1.5V stars. Nevertheless, 
this energy budget should be the lower limit because it is possible that some fraction of the 
ionizing photons would escape from the envelope of the surrounding neutral material. From 
our analysis, it is plausible that the HII region and the bright dust rim around [DBS2003] 45 
are excited by the central cluster. 
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5.2. Triggered Star Formation 



Observations indicate that the spatia lly separated s ubgroups in given OB asso ciations 
show systematic differences in stellar age telaauw 1964 ). Elmegreen fc Lada ( 1977 ) suggest 
that the formation of these subgroups is triggered sequentially by the expansion of HII re- 
gions. Other scen arios of triggered star formation on different scales have also been proposed 
(jElmegreenl 1 19981 . and references therein). As we know, in spontaneous star formation the 
fragmentation and collapse of molecular clouds is due to the loss of internal hydrostatic 
equilibrium between gravity and supporting forces. On the other hand, external dynamical 
processes dominate the compression of molecular clouds and initiate the following fragmen- 
tation and collapse in triggered star formation scenarios. Trigger ed star formation is believed 
to be ubiquitous across star forming regions and galactic disks (jElmegreenl 1 1 9981 ) . 



Among the proposed scenarios, the "collect and collapse" process at the peripheries of 
HII regions is particularly i nteresting because it is thought to be a mechanism by which 
very massive stars can form (IWhitworth et al.lll994j . , hereafter W94). In this scenario, star 
formation is triggered by the propagation of an ionization-shock front into molecular clouds. 
As the shock front drives into the clouds, neutral material is accumulated and compressed 
in a layer behind the front. Eventually, this layer of gas becomes gravitationally unstable, 
fragments, and collapses to form stars. Relatively higher temperature of the shocked gas 
causes the collapse to start at a higher Jean's mass. 

Despite its appealing to the theory of massive star formation, the collect and collapse 
process is not convincingly confirmed due to the lack of evidence. This is partially because 
young stellar objects in the peripheries of HII regions are usually embedded in molecular 
material and can not be easily detected at optical wavelengths. It is also difficult to dis- 
tinguish it from other triggering mechanisms since they all introduce systematic differences 
of the physical properties, e.g. stellar age and mass, between the formed subgroups. With 
appropriate initial conditions, even spontaneous star formation can cause similar observa- 
tional features. The initial hierarchical structure in large star forming regions can result in a 
systematic progression in size, age, and ve locity dispersion for OB subgroups formed through 
the spontaneous process (jElmegreenlll998l ) . Particularly, the collect and collapse process and 
the external pressure induced collapse of pre-existing clumps occur on the same scale and 
cause radial progressions in physical properties, therefore are difficult to be distinguished 
from each other. 



Elmegreen! (119981 ) points out that the collect and collapse process differs from the clump- 



squeezing process in two aspects. In the former scenario, star formation activities are more 
like a burst along the ridge of the accumulated gas. In contrast, star formation is continuous 
throughout the disturbed cloud in the latter case. The collect and collapse scenario tends 
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to form stellar clusters and no star at an intermediate age should exist between the clus- 
ters, while the clump compression mechanism tends to form individual stars or small stellar 
systems and the age distribution of stars is smooth. Therefore, it is possible to distinguish 
two processes using these properties. We should note that different mechanisms of star for- 
mation may be present at the same time. The shock front of an HII region can compress 
pre-existing globules inside molecular clouds and force them to collapse. At the same time, 
ambient medium between the globules is swept-up into a layer which can form stars through 
the collect and collapse process. Nevertheless, observations of young stars along the ridges 
of HII regions provide supporting evidence of star formation through the collect and collapse 
mechanism since stars formed through other mechanisms will not distribute preferentially 
along the shock front. 



In a search for infrared sources at HII region peripheries, iDeharveng et al.l (120051 ) ob- 
served bright and red objects in the dust rim surrounding the ionized gas around [DBS2003] 45. 
These authors suggest that the point sources are the results of the collect and collapse star 
formation. Our new observations of [DBS2003] 45 with Spitzer IRAC show a ring-like struc- 
ture in all four wavelength channels surrounding the cluster. Many bright sources can be 
seen along the ridge of the structure (Figure [TU| right). The light in these mid-IR channels 
i s generally attribute d to large molecules such as polycyclic aromatic hydrocarbons (PAHs) 



( ILeger &: Pugetlll984j ) . The observed morphology clearly indicates that the ring-like structure 
is the boundary of the ionized gas and very likely consists of neutral material accumulated 
during the expansion of the ionization front. Bright and red stellar objects spotted along 
this structure are the best example of triggered star formation in the swept-up layer around 
HII regions. Therefore, our observations support that the collect and collapse star formation 
process is undergoing in [DBS2003] 45. Further support for this star formation mechanism 
in action would come from observations of the age and mass distributions of the stars along 
the HII region peripheries. These stars should be significantly younger than the stars be- 
longing to or near the cluster. Observations of similar masses and separations for the newly 
formed stars along the peripheries would support that they are formed together through the 
fragmentation of the compressed layer since th is process occurs on a typical length scale and 



produces molecular cores with a typical mass (jWhitworth et al.lll994l ; lElmegreenl 119981 ) 



5.3. Comparison With A Theoretical Model 



Whitworth et al.l (119941 ) studied the fragmentation of the dense layer of accumulated 



material surrounding an expanding HII region and derived formula for the time at which the 
fragmentation starts (tf rag ), the radius of the HII region at that time (R/ ms ), the column 
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density of the layer, the masses of the fragments (Mf rag ) and the separation between them. 
Particularly, W94 showed that both tf rag and R/ ra9 depend weakly on the isothermal sound 
speed (&s) in the compressed layer, the number of Lyman continuum photons emitted by 
the exciting star per second (Q,Ly) an d the density of the ambient medium (n ): 



Rfrag 
tfrag 



5.8 a . 



4/11 n l/ll -6/11 



Q 



19 



1.56 a 7 / 11 Q]sa a. 



n 3 

-1/11 



pc 

-5/11 



'49 



Myr, 



(9) 
(10) 



where a.2=[as/0.2 km s _1 ], Q49=[Ql?//10 49 s _1 ] and n 3 =[n /10 3 cm -2 ] are dimensionless 
variables. They also predicted that the mass of fragments depends strongly on the sound 
speed ag: 



M 



frag — 23 CI 2 ^ Q49 ^ 



(11) 



This model has been tested by iDale et al.l ( 20071) usin g an 1-D smooth particle hydrody- 



namical code. Although Mf rag predicted by lDale et al.l is ~2.5 times smaller than the value 
predicted by the W94 model, the predicted tf rag and R/ ra9 from two models agree with each 
other within 20 to 25 per cent, suggesting that the W94 model is valid and accurate. 



Zavagno et al.l (120061 ) (hereafter Z06) applied the analytical mo del to a real HII region 



RCW 79 to estimate its age and other properties. Located at 4.2 kpc (jRusseilll2003l ). RCW 79 
appears to be a circular HII region surrounded by a partially open dust ring. Point sources 
are observed along the dust ring. This morphology of RCW 79 suggest that the expansion of 
the HII region tr iggered the star format ion in the dust ring as described by the collect and 
collapse process (IDeharveng et al.ll2005l ). Z06 compared the 1.2 mm continuum, NIR and 
mid-IR data of RCW 79 with the analytical model and estimated that the dynamical age of 
RCW 79 is 1.7 Myr. They also estimated that the triggered star formation along the dust 
ring occurred ~10 5 years ago. These timescales and the corresponding density (~2000 cm -3 ) 
of the surrounding molecular material are used to constrain the age of an associated compact 
HII region and the masses of fragments observed in the ring. The good agreement between 
the observations and the model predictions indicates that the collect and collapse process is 
in deed the main triggering mechanism of star formation in the dust ring around RCW 79. 

We have not performed molecular line observations to help us constrain the masses of 
fragments and the sound speed inside the molecular layer in [DBS2003] 45. However, using 
the morphology of [DBS2003] 45 and the theoretical predictions of W94 (similar to the study 
of RCW79 by Z06), we can compare the fragmentation age of the expanding bubble with both 
the bubble's dynamical age and the age derived for the cluster. The Spitzer IRAC images 
of [DBS2003] 45 show that the average radius of the ring-like structure is about 8 arcmin. 
This corresponds to a fragmentation radius of <10.5 pc at a distance of 4.5 kpc. Adopting 
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an intermediate value (0.4 km s -1 ) for as, we can derive no ~500 cm -3 and tf rag ~3.3 Myr, 
based on Eq. [9] and [10J Since we assume that the molecular cloud is uniform with a small 
sound speed and the expansion of the HII region is continuous, the derived n should be 
considered as the lower limit. Here we have applied Q^,,=8 .2xl0 48 photons per second. 



Additionally, according to W94 and iDyson fc Williams! (119971 ). the radius of the swept-up 



shell around an evolved HII region grows with time t: 

fl~4.5 QH 7 n~ 2/7 t% 7 pc, (12) 

where tM=[t/Myr]. Adopting R=10.5 pc, Q 49 =0.82 and n 3 =0.5, we have t ~3.3 Myr. The 
above age estimates are consistent with each other and suggest that the point sources along 
the dust ridge formed not long time ago. A typical sound speed in the compressed layer 
is around 0.2-0.6 km s _1 . We used as=0.4 km s _1 in the above calculation, which implies 
an uncertainty of ±1 Myr for the derived fragmentation age. We have also assumed other 
values, e.g. the density of the compressed layer, the ionizing luminosity of the central cluster 
and the size of the ionized bubble, which may introduce additional uncertainty into the 
derived dynamical age and the fragmentation time of the bubble. If we take the uncertainty 
into account, 3.3 Myr is in a good agreement with the age of the central cluster of 5 Myr 
that we have derived in the earlier section. 

The above derived numbers rely on the following assumptions: the ionizing photon 
luminosity of the exciting source is constant during the lifetime of the HII region; the HII 
region expands within an originally uniform and smooth molecular cloud; the compressed 
layer by the ionization/shock front has a typical sound speed in dense molecular clouds; 
the effects of stellar winds and other more complex physical processes are not considered. 
These assumptions may be oversimplified compared to real star forming regions. More 
rigorous calculation of the nebular dynamics is well beyond the scope of the current work. 
Nevertheless, our simple calculations suggest that the radio emission from the HII region is 
consistent with being powered by the ionizing flux from the cluster. The dynamical age and 
the fragmentation age of the expanding bubble are consistent with the age we derived for the 
central cluster. These findings together point to a scenario in which the central cluster has 
driven the expansion of the HII region into the surrounding medium, and a second generation 
of star-formation is currently being triggered through the collect and collapse mechanism at 
the periphery of the bubble. 



6. Summary 



We present a photometric and spectroscopic study of a stellar cluster candidate on the 
southern sky. Our study confirms that [DBS2003] 45 is a cluster of a total mass ~10 3 M Q . 
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The field decontaminated CMD clearly shows a cluster sequence at H-Ks~0.45, suggesting 
an average extinction Av~7.1 Mags. Based on low-resolution spectroscopy, we identify a few 
bright cluster members to be early type blue giants/supergiants. The presence of these late 
O / early B type stars and the absence of any early O type stars and M-type supergiants, 
suggests an age between 5 and 8 Myr for the identified cluster. A mass function with an 
approximately Salpeter slope is found for the cluster. The associations of the cluster with the 
SuperCOSMOS H a emission morphology and the ring-shaped 8.0/iin structure observed by 
the infrared camera IRAC on the Spitzer telescope support that the cluster is relatively young 
and the ambient medium has not been completely dispersed by stellar winds and supernova 
explosions. Evidence of patchy extinction toward the cluster is found and is consistent with 
the presence of a population of heavily reddened stars. The presence of infrared point sources 
along the rim of the dusty bubble supports the mechanism of triggered star formation. Simple 
calculations of the radio emission and the ages of the HII region support that the identified 
cluster is the source to drive the expansion of the HII region into the ambient medium and 
to trigger the formation of new stars. 

The material in this paper is supported by NASA under award NNG 05-GC37G, through 
the Long Term Space Astrophysics program. This research was performed in the Rochester 
Imaging Detector Laboratory with support from a NYSTAR Faculty Development Program 
grant. The work is also supported by a special grant from National Natural Science Foun- 
dation of China, NSFC- 10843008. Part of the data presented here was obtained with the 
ESO New Technology Telescope (PI: Benjamin Davies, PID: 080.D-0470). This publication 
makes use of data products from the Two Micron All Sky Survey, which is a joint project of 
the University of Massachusetts and the Infrared Processing and Analysis Center, Califor- 
nia Institute of Technology, funded by the National Aeronautics and Space Administration 
and the NSF. This research has made use of Spitzer's IRAC imaging data (PI: Donald 
Figer, PID:30734), superCOSMOS's H a line sky survey and Spitzer's massive star cluster 
servery and galactic structure and star-forming region imaging data (PI: Steven R Majewski, 
PID:40791), PMN radio 3842 MHz sky survery data, the SIMBAD database, and the GSFC 
IDL library. We also thank Dr. Lise Deharveng for advices on radio observations. 

REFERENCES 

Bica, E., Dutra, C. M., Soares, J., & Barbuy, B. 2003, A&A, 404, 223 
Blaauw, A. 1964, ARA&A, 2, 213 

Blum, R. D., Conti, P. S., & Damineli, A. 2000, AJ, 119, 1860 



-19- 



Dale, J. E., Bonnell, I. A., & Whitworth, A. P. 2007, MNRAS, 375, 1291 

Deharveng, L., Zavagno, A., & Caplan, J. 2005, A&A, 433, 565 

Dias, W. S., Alessi, B. S., Moitinho, A., & Lepine, J. R. D. 2002, A&A, 389, 871 

Dutra, C. M., Bica, E., Soares, J., & Barbuy, B. 2003a, A&A, 400, 533 

Dutra, C. M., Ortolani, S., Bica, E., Barbuy, B., Zoccali, M., & Momany, Y. 2003b, A&A, 
408, 127 

Dyson, J. E. & Williams, D. A. 1997, The physics of the interstellar medium, ed. J. E. Dyson 
& D. A. Williams 

Elmegreen, B. G. 1998, in Astronomical Society of the Pacific Conference Series, Vol. 148, 
Origins, ed. C. E. Woodward, J. M. Shull, & H. A. Thronson, Jr., 150-+ 

Elmegreen, B. G. & Lada, C. J. 1977, ApJ, 214, 725 

Figer, D. F., Kim, S. S., Morris, M., Serabyn, E., Rich, R. M., & McLean, I. S. 1999, ApJ, 
525, 750 

Froebrich, D., Scholz, A., & Raftery, C. L. 2007, MNRAS, 374, 399 

Hanson, M. M., Kudritzki, R.-P, Kenworthy, M. A., Puis, J., & Tokunaga, A. T. 2005, 
ApJS, 161, 154 

Kim, S. S., Figer, D. F., Kudritzki, R. P., & Najarro, F. 2006, ApJ, 653, L113 

Kroupa, P. 2001, MNRAS, 322, 231 

-. 2002, Science, 295, 82 

Leger, A. & Puget, J. L. 1984, A&A, 137, L5 

Lejeune, T. & Schaerer, D. 2001, A&A, 366, 538 

Moorwood, A., Cuby, J.-G., & Lidman, C. 1998, The Messenger, 91, 9 
Morris, M. 1993, ApJ, 408, 496 

Parker, Q. A., Phillipps, S., Pierce, M. J., Hartley, M., Hambly, N. C, Read, M. A., 
MacGillivray, H. T., Tritton, S. B., Cass, C. P., Cannon, R. D., Cohen, M., Drew, 
J. E., Frew, D. J., Hopewell, E., Mader, S., Malin, D. F., Masheder, M. R. W., 
Morgan, D. H., Morris, R. A. H., Russeil, D., Russell, K. S., & Walker, R. N. F. 2005, 
MNRAS, 362, 689 



-20- 



Rieke, G. H. & Lebofsky, M. J. 1985, ApJ, 288, 618 
Russeil, D. 2003, A&A, 397, 133 

Rybicki, G. B. & Lightman, A. P. 1986, Radiative Processes in Astrophysics (Radiative 
Processes in Astrophysics, by George B. Rybicki, Alan P. Lightman, pp. 400. ISBN 
0-471-82759-2. Wiley- VCH , June 1986.) 

Salpeter, E. E. 1955, ApJ, 121, 161 

Scalo, J. 1998, in Astronomical Society of the Pacific Conference Series, Vol. 142, The Stellar 
Initial Mass Function (38th Herstmonceux Conference), ed. G. Gilmore & D. Howell, 
201-+ 

Schaerer, D. & de Koter, A. 1997, A&A, 322, 598 

Schaller, G., Schaerer, D., Meynet, G., & Maeder, A. 1992, A&AS, 96, 269 

Sternberg, A., Hoffmann, T. L., & Pauldrach, A. W. A. 2003, ApJ, 599, 1333 

Stolte, A., Grebel, E. K., Brandner, W., & Figer, D. F. 2002, A&A, 394, 459 

Whitworth, A. P., Bhattal, A. S., Chapman, S. J., Disney, M. J., & Turner, J. A. 1994, 
MNRAS, 268, 291 

Wright, A. E., Griffith, M. R., Burke, B. F., & Ekers, R. D. 1994, ApJS, 91, 111 

Zavagno, A., Deharveng, L., Comeron, F., Brand, J., Massi, F., Caplan, J., & Russeil, D. 
2006, A&A, 446, 171 

Zinnecker, H. & Yorke, H. W. 2007, ARA&A, 45, 481 



This preprint was prepared with the AAS IATgX macros v5.2. 



- 21 - 




Fig. I. — Ks band mosaic of the cluster candidate [DBS2003] 45 in a logarithmic scale. The 
origin of the axes corresponds to R.A.=10 /l 19 m 10.5 s , Decl.=-58°02'22.6" (2000). 
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Fig. 2. — Using 2MASS magnitudes to calibrate the observed photometry. In both figures, 
plus signs indicate H band magnitudes and diamonds in the figures indicate Ks band mag- 
nitudes. In the left figure, x and y axes are observed and 2MASS magnitudes, respectively. 
Solid lines and equations show the results of a single variable linear model fitting. In the right 
figure, x axis is the calibrated magnitude and y axis is the difference between the calibrated 
and the instrumental magnitudes. 
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Fig. 3. — The completeness curves of the derived photometry for the candidate cluster when 
choosing different sizes. The curve for the outermost region shows that the photometry is 
50% complete at Ks~17. The dashed line indicates 50% completeness. 
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Fig. 4. — Low-resolution spectra of a few picked stars. The image shows a blowup of the 
Ks band mosaic of the cluster with symbols indicating different groups of stars in the NIR 
CMD (Figure [7]) and the mid-IR color-color diagram (Figure [9]) of the cluster. Filled circles 
represent stars detected by SofI and stay in the CMD after contamination subtraction. 
Squares represent stars detected by IRAC. Diamonds and numbers indicate locations of 
stars with a spectrum shown in the right panel. The image is centered at (RA, DEC) = 
(10 h 19 m 10.5 s , -58°02'22.6"). 
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Fig. 5.— The Color-Magnitude diagram (CMD) of [DBS2003] 45. The Ks band image 
(in a logarithmic scale) on the left shows the region from which the CMDs are drawn. 
Two concentric circles centered at (RA, DEC) = (10H9 m 10.5 s , -58°02'22.6") demonstrate 
the central circular region and an annulus region chosen for the cluster candidate and the 
control field. The radii of the circles are 0.72 and 2.4 arcmin, respectively. The CMDs on 
the right in the figure are for the cluster with background contamination (left panel), the 
decontaminated cluster (middle panel) and the control field CMD (right panel), respectively. 
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Fig. 6. — The number of stars with masses over 15 M versus the age in a cluster of mass 10 3 , 
2xl0 3 and 10 4 M Q , respectively. The 'Salpeter' IMF is assumed for the mass distribution 
at the zero age. The Geneva isochrones with solar abundances are used for the evolution of 
cluster members. The solid straight line indicates the number of massive stars detected in 
[DBS2003] 45. 
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Fig. 7. — The CMD of [DBS2003] 45 with field contamination removed. The solid line shows 
the 6 Myr isochrone. The stars with low resolution spectra are signified with diamonds 
and numbers. The locations of stars of certain initial masses are also indicated with the 
reddening vector. 
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Fig. 8. — (Left) The mass function of [DBS2003] 45 assuming a 6 Myr isochrone. The 
dashed histogram shows the result before completeness correction and the solid histogram 
shows the result after the correction. The fitting result to a power law is indicated for the 
solid histogram. A classic 'Salpeter' model is also plotted. (Right) The slope of the mass 
function with different bin numbers. 




Fig. 9. — Left: The IRAC color-color diagrams of the cluster. Squares represent data points 
in the IRAC photometry and filled circles represent the stars with SofI measurements. Stars 
with spectra are indicated with diamonds and numbers. Right: The locations of stars 
showing 5.8 fxm excess (Group II) and strong foreground dust extinction (Group III) are 
indicated with squares and diamonds, respectively. The image is centered at (RA, DEC) = 
(10 h 19 m 10.5 s , -58°02'22.6"). 
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Fig. 10.— Left: The contours of the IRAC 8.0 fim image are plotted over the SUPERCOS- 
MOS Hq, grey scale image. It is evident that emission from neutral material is mainly come 
from a ring-like structure around the concentration of emission of ionized gas identified by 
H a line emission. The location of the cluster is indicated with a circle. Right: The IRAC 
8.0 jum image is shown in grey scale. Point sources can be seen at several locations along the 
rim of emission. Both images are centered at (RA, DEC) = (10H9 m 10.5 s , -58 o 02'22.6"). 
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Table 1. Magnitude Calibration Coefficients 



Band 


a(A) 


b(B) 


c(C) 


h 


2.396 


1.008 


0.007 


ks 


2.876 


0.001 


1.019 


H 


-2.359 


0.992 


-0.007 


Ks 


-2.822 


-0.001 


0.982 



Note. — Magnitude transforma- 
tion coefficients derived through the 
regression between 2MASS magni- 
tudes and instrumental magnitudes. 
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Table 2. LIST OF STARS 



Star# RA(J2000) DEC(J2000) 



1 


154 


48 


00.5 


- 58 


02 


23.4 


2 


154 


47 


31.7 


- 58 


02 


26.6 


3 


154 


47 


15.6 


- 58 


02 


34.6 


4 


154 


47 


47.8 


- 58 


02 


13.0 


5 


154 


47 


13.8 


- 58 


02 


13.9 


6 


154 


47 


35.7 


- 58 


01 


59.5 


7 


154 


17 


26.5 


- 58 


01 


55.5 


8 


154 


46 


39.9 


- 58 


01 


16.0 



Note. — Coordinates of the stars indicated 
in Figure [4] 



